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Abstract 

The 3-D supersonic viscous flow in an experimental 
MHD channel has been numerically simulated. The 
experimental MHD channel is currently in operation 
at NASA Ames Research Center. The channel con- 
tains a nozzle section, a center section, and an ac- 
celerator section where magnetic and electric fields 
can be imposed on the flow. In recent tests, veloc- 
ity increases of up to 40% have been achieved in the 
accelerator section. The flow in the channel is numer- 
ically computed using a new 3-D parabolised Navier- 
Stokes (PNS) algorithm that has been developed to 
efficiently compute MHD flows in the low magnetic 
Reynolds number regime. The MHD effects are mod- 
eled by introducing source terms into the PNS equa- 
tions which can then be solved in a very efficient man- 
ner. To account for upstream (elliptic) effects, the 
flowfield can be computed using multiple streamwise 
sweeps with an iterated PNS algorithm. The new al- 
gorithm has been used to compute two test cases that 
match the experimental conditions. In both cases, 
magnetic and electric fields are applied to the flow. 
The computed results are in good agreement with the 
available experimental data. 

Introduction 

Magnetohyrodynamics (MHD) can be utilized to im- 
prove performance and extend the operational range 
of many systems. Potential applications include hy- 
personic cruise, advanced Earth-to-orbit propulsion, 
chemical and nuclear space propulsion, regenerative 
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aerobraking, onboard flow control systems, test fa- 
cilities, launch assist, and power generation. One of 
the critical technologies associated with these appli- 
cations is MHD acceleration. In order to study MHD 
acceleration, an experimental MHD channel has been 
built at NASA Ames Research Center by D. W. Bog- 
danoff, C. Park, and U. B. Mehta [1,2]. The channel 
is about a half meter long and contains a nozzle sec- 
tion, a caiter section, and an accelerator section. The 
channel has a uniform width of 2.03 cm. Magnetic 
and electric fields can be imposed upon the flow in 
the accelerator section. A cross section of the MHD 
ch ann el is shown in Fig. 1. 

In the present study, the flow in the experimen- 
tal MHD channel is numerically simulated. Flow- 
fields involving MHD effects have typically been 
computed [3-15] by solving the complete Navier- 
Stokes (N-S) equations for fluid flow in conjunction 
with Maxwell’s equations of electromagnetodynam- 
ics. When chemistry and turbulence effects are also 
included, the computational effort required to solve 
the resulting coupled system of partial differential 
equations is extremely formidable. One possible rem- 
edy to this problem is to use the parabolized Navier- 
Stokes (PNS) equations in place of the N-S equar- 
tions. The PNS equations can be used to compute 
three-dimensional, supersonic viscous flowfields in a 
very efficient manner [16]. This efficiency is achieved 
because the equations can be solved using a space- 
marching technique as opposed to the time-marching 
technique that is normally employed for the complete 
N-S equations. 

Recently, the present authors have developed PNS 
codes to solve 2-D supersonic MHD flowfields in both 
the high and low magnetic Reynolds number regimes 
[17,18]. The magnetic Reynolds number is defined 
as Rc-rn — <7 e p c Foo£ where a c is the electrical con- 
ductivity, fj> c is the magnetic permeability, Fqo is the 
freestream velocity, and L is the reference length. 
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The new MUD PNS codes are based on NASA’s up- 
wind PNS (UPS) code which was originally devel- 
oped by Lawrence et al. [19]- The UPS code solves 
the PNS equations using a fully conservative, finite- 
volume approach in a general nonorthogonal coordi- 
nate system. The UPS code has previously been ex- 
tended to permit the computation of flowfields with 
strong upstream influences- In regions where strong 
upstream influences are present, the governing equa- 
tions are solved using multiple sweeps (i.e. itera- 
tions). As a result of this approach, a complete flow- 
field can be computed more efficiently (in terms of 
computer time and storage) than with a standard N- 
S solver which marches the entire solution in time. 
Three iterative PNS algorithms called IPNS, TIPNS, 
and FBIPNS have been developed and are described 
in Refe. [20-22]. 

For many aerospace applications, including the 
present experimental MHD channel, the electrical 
conductivity of the fluid is low and hence the mag- 
netic Reynolds number is small. In these cases, it 
makes sense to use the low magnetic Reynolds num- 
ber assumption and reduce the complexity of the gov- 
erning equations. The MHD effects are modeled with 
the introduction of source terms into the fluid flow 
equations* as was done in the present low magnetic 
Reynolds number PNS code [18]. This code has been 
extended to three-dimensions in the present study 
and is used for the numerical simulations of the flow 
in the experimental MHD channel. 


Governing Equations 

The governing equations for a viscous MHD flow 
with a small magnetic Reynolds number are given 
by [14]: 

Continuity equation 

|? + V(pV) = 0 (1) 

Momentum equation 

^+V-[pW+pi]=V-f + JxB (2) 
Energy equation 

^p+V-[(pe* +1>) V] = V-(V-f)— V-U+B-J (3) 
Ohm’s law 

J = cr e (E + VxB) (4) 

where V is the velocity vector, B is the magnetic 
field vector, E is the electric field vector, and J is the 
conduction current density. 


The governing equations are nondimensi onalized 
using the following reference variables: 
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where the superscript * refers to the nondimen- 
sional quantities. For convenience, the asterisks are 
dropped in the following equations. 

The governing equations written in vector form in 
a 3-D Cartesian coordinate system become 


OTJdEidTidGi _ dE v dF„dG v i c 
dt + dx + dy + dz dx + dy + dz +Smhd 

( 6 ) 

where U is the vector of dependent variables, E*, F » 
and G{ are the in viscid flux vectors, and E w , F„ and 
G v are the viscous flux vectors. The source term 
Smhd contains all of the MHD effects. The flux vec- 
tors are given by 


U = [ p, pu, pv, pw, pet ] T (7) 
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(11) 
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+E Z {E Z + u By — vB z ) 


where 


pet = ip^ + t^ + ur 2 ) + 


(12) 


(13) 


and the nondimensional shear stresses and heat fluxes 
are defined in the usual manner [16]- 

The flow can be computed assuming either a con- 
stant 7 or by using the simplified curve fits of Srini- 
vasan et al. [23,24] for the thermodynamic and trans- 
port properties of equilibrium air. For the latter case, 
7 in Equation (13) is determined using the curve fit 
7 = 7(e,p). In future calculations, the flow will be 
computed in chemical nonequilibrium. 

The governing equations are transformed into com- 
putational space and written in a generalized coordi- 
nate system (f , rj, C) as 


iu t + E e +F„ + G c = ^p (14) 

where 

E = (y) (E< - E„) + (&) (f 4 - F.) 

+ ( 7 ) (g ‘- g ’> 

F = (jj) (E» — E«) + (jj) (Fi — F„) 

+ { T j) (G<_G " ) 

G = ^(E,-E.)+0^(F 4 -F.) 

+ ^y^ ( G « ~ Gw ) (15) 

and J is the Jacobian of the transformation. 

The governing equations are parabolized by drop- 
ping the time derivative term and the streamwise 


direction (f) viscous flow terms in the flux vectors. 
Equation (14) can then be rewritten as 

E ? + F„ + G c = ^p (16) 

where 

* - (tt) ^ + (!f) Fi + ( /) G< 

F = (^)(B 1 -E'J + (i)(F < -r v ) 

+ (t) (G<_g;) 

g = ( Ei -E;) + ^) (Fi-r v ) 

+ (y) (G < - G t) (17) 

The primes in the preceding equations indicate that 
the streamwise viscous flow terms have been dropped. 

For turbulent flows, the two-layer Baldwin-Lomax 
turbulence model [25] has been modified to account 
for MHD effects. Only the expression for turbulent 
viscosity in the inner layer is changed. This modifi- 
cation for MHD flows is due to Lykoudis [26]. 


Numerical Method 

The governing PNS equations with MHD source 
terms have been incorporated into NASA’s upwind 
PNS (UPS) code [19]. These equations can be solved 
very efficiently using a single sweep of the flowfield 
for many applications. For cases where upstream 
(elliptic) effects are important, the flowfield can be 
computed using multiple streamwise sweeps with ei- 
ther the IPNS [20], TIPNS [21], or FBIPNS [22] algo- 
rithms. This iterative process is continued until the 
solution is converged. 

For the iterative PNS (IPNS) method, the E vec- 
tor is split using the Vigneron parameter (uj) [27]. 
This parameter does not need to be changed for the 
present low magnetic Reynolds number formulation. 
In the previous high magnetic Reynolds number code 
[17] it was necessary to modify the Vigneron param- 
eter to account for MHD effects. After splitting, the 
E vector can be written as: 

E = E*+E p (18) 

where 
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The streamwise derivative of E is then differenced 
using a backward difference for E* and a forward dif- 
ference for the “elliptic” portion (E J> ): 


(^) = ^ [(®<+i - E*) + (Ef +2 - Ef +1 )] 

( 20 ) 

where the subscript (t 4- 1) denotes the spatial index 
(in the f direction) where the solution is currently 
being computed. The vectors E* +1 and E^ are then 
linearized in the following manner: 


E; +1 = E*+^^(U, +1 -U t ) 

Ef+i = Ef+(^) j (U i +i-U t ) 
The Jacobians can be represented by 


( 21 ) 
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( 22 ) 


After substituting the above linearizations into 
Eq. (20), the expression for the streamwise gradient 
of E becomes 


+ (Ef+2 ~ Ef) ] (23) 


The final discretized form of the fluid flow equa- 
tions with MHD source terms is obtained by substi- 
tuting Eq. (23) into Eq. (16) along with the linearized 
expressions for the fluxes in the cross flow plane. The 
final expression becomes: 


±(a--a?)+ 9 < dF 


.d_(9G 

• r\ > 


-I fc+1 


dv \dv j, ac 

= RHS 


(AUi)* +1 

(24) 


where 


(AUi)* +1 = (U i+1 -U i )‘ +1 



and the superscript k+1 denotes the current itera- 
tion (i.e. sweep) level. In the preceding equation, the 
MHD source term, Smhd, is treated explicitly since it 
is evaluated using the velocity at station i (V*). For 
most cases, this will not degrade the accuracy of the 
solution since Af is small and the velocity changes 
slowly. If this is not the case, a predictor-corrector 
procedure can be implemented whereby a predicted 
velocity at station i+1 (V\~ ) is first obtained using 
Eq. (24). The solution at station i-Hl is then recom- 
puted by evaluating Smhd with V~ - 


Numerical Results 

The numerical calculation of the supersonic flow in 
the experimental MHD channel is now discussed. 
This flowfield was previously computed by MacCor- 
mack [12] using the full N-S equations coupled with 
the electromagnetodynamic equations. The electrical 
conductivity in his calculations was set at 1.0 x 10 4 5 
mho/m resulting in a very large magnetic Reynolds 
number. In the present study, the calculations are 
performed in the low magnetic Reynolds number 
regime using a realistic value of electrical conductiv- 
ity. Both 2-D and 3-D results have been obtained. 

The flow in the nozzle section was computed using 
a combination of the OVERFLOW code [28] and the 
present PNS code (without MHD effects). For the 3- 
D OVERFLOW nozzle calculation, a highly stretched 
grid consisting of 130 x 50 x 50 grid points was used. 
The normal grid spacing at the wall was 1.0 x 10“ 5 
m. For the PNS calculation of the flow in the re- 
mainder of the nozzle and the rest of the MHD chan- 
nel, a highly stretched grid consisting of 60 points in 
both the y and z directions was used and the nor- 
mal grid spacing at the wall was 1.0 x 10“ 5 m. As 
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a consequence of flow symmetry, only one-fourth of 
the channel cross section was computed in the 3-D 
calculations. 

The calculations were performed assuming turbu- 
lent flow throughout the MHD channel. In addition 
the flow was computed using a constant 7 of 1.25 
to simulate equilibrium air. The channel wall tem- 
perature was assumed to be isothermal since steady 
flow conditions were maintained in the experiment for 
only about 1.2 milliseconds. A schematic of the pow- 
ered portion of the MHD channel with the directions 
of the applied magnetic and electric fields is shown 
in Fig. 2. The values of the electrical conductivity 
(cr e ), the magnetic field (jB z ), and the electric field 
(E y ) were kept constant in the powered portion of 
the channel. Two test cases corresponding to Runs 
15 and 16 of the NASA Ames experiments [29] were 
computed in this study and are now discussed. 

Test Case 1: NASA Ames MHD Run 15 

(Vcap. = 32QV) 

The dimensional flow parameters for this test case 
are: 

Po = 9.10 x IQ 5 N/m 2 

T 0 = 5560 K 

T w = 300 K 

<r e ~ 130 mho/m 

B z = 0.0, 0.92 T 

E y = 0, 3955, 5000, 6000 Y/m 

where the subscript o denotes total conditions at the 
nozzle entrance and w denotes wall conditions. 

This case was computed using several different elec- 
tric field strengths in order to properly simulate the 
experiment. In the experiment, the voltage applied 
to the electrodes was approximately 134 V for this 
case, however, due to the sheath voltage drop, the 
actual voltage applied to the flow is smaller than the 
electrode voltage. The voltage drop was measured for 
the central inviscid core flow, and was approximately 
67 V [2]. Since the boundary layer is computed in the 
numerical solution, the applied electric field must be 
approximately the voltage drop across the electrodes 
minus the sheath voltage drop. Unfortunately, it is 
not a trivial task to measure the sheath voltage drop. 
Therefore, several different electric fields were cho- 
sen in the numerical calculations so that the corre- 
sponding voltage drop across the electrodes would be 
between 67 V and 134 V. The voltage drop of 67 V 
corresponds to E y = 3955 V/m at the center of the 
accelerator section and a voltage drop of 101.6 V cor- 
responds to E y = 6000 V/m, with a sheath voltage 
drop of 32.4 V. 


The computed streamwise variation of static pres- 
sure for the 2-D calculations is shown in Fig. 3 for the 
different electric field strengths. The pressure varia- 
tion with no electric field or magnetic field is denoted 
by E y = 0. The results are in reasonable agreement 
with the experiment. The numerical results show an 
increase in static pressure as the electric field strength 
is increased. 

The computed streamwise variation of averaged ve- 
locity for the 2-D calculations is shown in Fig. 4. The 
velocities are averaged across the channel cross sec- 
tion and normalized using the entrance velocity to be 
consistent with the experiment. In the experiment, 
the velocities were obtained by measuring the volt- 
age generated by the flow at the last electrode pair 
(19) which is unpowered. This procedure inherently 
involves an averaging of the velocity profile. The nu- 
merical results indicate an increase in the averaged 
velocity of about 26% with E y = 6000 V/m and this 
compares well with the experimental value of approx- 
imately 27%. The velocity vector plots for Ey = 0 
and E y — 6000 V /m are shown in Fig. 5. 

The computed streamwise variation of averaged ve- 
locity for the 3-D calculations is shown in Fig. 6. The 
velocities are averaged and normalized in the same 
manner as for the 2-D calculations- The numerical 
results indicate an increase in the averaged velocity 
of about 26% with E y = 6000 V/m which is the same 
value obtained with the 2-D calculations. 

Test Case 2: NASA Ames MHD Run 16 

(V cap . = 380V) 

The dimensional flow parameters for this test case 
are: 

Po = 9.92 x10 s N/m 2 
To = 5560 K 
T w — 300 K 
a e = 140 mho/m 
B z = 0.0, 0.92 T 

E y = 0, 4309, 5000, 6000, 7000 V/m 

This test case was also computed using several dif- 
ferent electric field strengths in order to properly 
simulate the experiment. The computed stream- 
wise variation of static pressure for the 2-D calcula- 
tions is shown in Fig. 7 for the different electric field 
strengths. The computed pressures are in good agree- 
ment with the experimental pressures. The computed 
streamwise variation of averaged velocity for the 2-D 
calculations is shown in Fig. 8. The numerical results 
indicate an increase in the averaged velocity of about 
31% with E y = 7000 V /m. This is less than the value 
of 38% that was obtained in the experiment. The cor- 
responding streamwise variation of averaged velocity 
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for the 3-D calculations is shown in Fig. 9. The nu- 
merical results indicate an increase in the averaged 
velocity of about 30% with E y = 7000 V /m which is 
again less than the experimental value of 38%. 

Concluding Remarks 

In this study, a new 3-D parabolized Navier-Stokes 
algorithm has been developed to efficiently compute 
MHD flows in the low magnetic Reynolds number 
regime. The new algorithm has been used to com- 
pute the flow in the NASA Ames experimental MHD 
channel for Runs 15 and 16. The numerical results 
are in good agreement with most of the experimental 
results. 
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Figure 2: Schematic of powered portion of MHD channel 
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Figure 4: Streamwise variation of averaged velocity (2-D) 
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Figure 6: Streamwise variation of averaged velocity (3-D) 


10 

American Institute op Aeronautics and Astronautics 
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Figure 7: Streamwise variation of static pressure (2-D) 



Figure 8: Streamwise variation of averaged velocity (2-D) 
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Figure 9: Stream wise variation of averaged velocity (3-D) 


12 

American Institute op Aeronautics and Astronautics 



